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METHOD AND APPARATUS FOR BACKSCATTER SPECTROSCOPY 

The present invention relates to a method of determining of a physical feature of a 
5 medium, comprising: 

- producing radiation with a light source; 

- placing a probe on a sample of said medium, the probe comprising a first optical fiber 
havuig a first diameter, and at least a second optical fiber having a second diameter; 

- sending light coming from the light source, through the first optical fiber; 

10 - collecting first backscattered radiation through the first optical fiber and second 
backscattered radiation through the second optical fiber; 

- producing a first signal based on the fibrst backscattered radiation, and a second signal 
based on the second backscattered radiation; 

- determine a measured differential backscatter signal as a function of wavelength using 
15 the first and second signals. 

Such a method is known firom Amelink et al [1]. There, a special device is used to 
determine particle sizes in superficial layers. The device is suitable for measuring 
particle sizes in for escample an aqueous suspension with polystyrene spheres, but is not 
fitted to accurately measure particle sizes in Uving tissue. So, determining whether 
20 living tissue is normal or precancerous, by way of measuring particle sizes hi hving 
tissue is not very promising. 

In Doombos et al [2] the optical properties of human tissue are determined in 
vivo using a spectroscopic arrangement with ten optical fibers. One of the fibers is used 
to irradiate a sample, and nine other fibers collect the reflected light By using a 
25 multitude of fibers to collect the reflected light, it is possible to calculate scattering and 
absorption coefficients of the sample. However, the method is not suitable for locally 
measuring the optical properties of the tissue. In particular, only mean values of the 
absorption coefficient of a relatively large part of the sample can be determined. 

It is an object of the present invention to locally measure a physical feature, such 
30 as a concentration, of a substance in a medium. 

The object is achieved by a method as described above, characterized by 

- calculating the physical feature by curve fitting said measured differential backscatter 
signal to a backscatter function, in which the backscatter fimction is a function of an 
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average pafh-length travelled by detected scattered photons, wherein the average path- 
length is independent from an absorption coefiBcient of the mediiim, and fiom a 
scattering coefficient of the medium. Contrary to methods using diffiisely scattered 
photons such as Doombos et al [2], in the method according to the invention, the local 

5 absorption coefficient of the sample is measured in an absolute way, independent of the 
magnitude of the local scattering and absorption coefficients. This facilitates the 
measurement of absolute concentrations of absorbing molecules in a sample without 
requiring prior knowledge of the magnitude of the scattering and absorption 
coefficients of the medium. 

10 In an embodiment, the average path-length is proportional to the first fiber 

diameter. This has as additional advantage that the average path-length and thereby the 
average penetration depth into the sample of the photons that contribute to the 
differential backscatter signal can be controlled by choosing the fiber diameter. As a 
result, the sampling volume can be controlled by adjusting the fiber diameter. 

1 5 Hence, the fiberoptic probe can be engineered to match the relevant dimensions of the 
medium under investigation. 

In a particular embodiment, the physical feature is a concentration of at least one 

substance in the medium. 

The invention also relates to a device for determining a physical feature of a 
20 medium, comprising: 

- a light source for producing radiation; 

- a probe with at least a first and a second optical fiber, the first optical fiber having a 
first diameter and being arranged to deliver the radiation on a sample of said medium 
and to collect first backscattered radiation fix)m said sample, the second optical fiber 

25 having a second diameter and being arranged to collect second backscattered radiation, 
wherein the second optical fiber is positioned alongside the first optical fiber; 

- a spectrometer for producing a first signal based on the first backscattered radiation, 
and for producing a second signal based on the second backscattered radiation; 

- a processor arranged to determine a measured differential backscatter signal as a 
30 fiinction of waveler^ using the first and second signals, 

characterized in that the processor is arranged to calculate the physical feature by curve 
fitting the measured differential backscatter signal to a backscatter fianction, hi which 
the backscatter fiinction is a fiinction of an average path-length travelled by detected 
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scattered photons, the average path-length being independent from an absorption 
coefficient of the medium, and from a scattering coefficient of the medium. 

Furthermore, the invention relates to a computer program according to claim 
Sand to a data carrier according to claim 9. 

5 

In another aspect of the invention, the invention relates to a method of 
determining a physical feature of a medium, comprising: 

- producing radiation with a Ught source; 

- placing a probe on a sample of the medium, the probe comprising a first optical fiber 
10 having a first diameter, and at least a second optical fiber having a second diameter; 

- sending light commg jfrom the light source, through the first optical fiber; 

- collecting first backscattered radiation through the first optical fiber and second 
backscattered radiation through the second optical fiber; 

- producing a first signal based on the first backscattered radiation, and a second signal 
1 5 based on the second backscattered radiation; 

- determine a measured differential backscatter signal as a function of wavelength using 
the first and second signals, 

characterized by 

- calculating the physical feature by curve fittmg the measured differential backscatter 
20 signal to a backscatter function, in which the backscatter function is a function of a 

mean firee path of photons. In this method, it is assumed that only singly scattered 
photons contribute to the differential backscatter signal and as a result the backscatter 
function can be easily derived analytically. 

In an embodiment, the physical feature is a concentration of at least one 
25 substance in the medium. 

The invention also relates to a device according to claim 13. 

Furthermore, the invention relates to a computer program according to claim 14 
and to a data carrier according to claim 15. 

Finally, the invention relates to a method according to claim 16. 

30 

The present invention will be described below with reference to exemplary 
embodiments and the accompanying schematic drawings, 'm which: 
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Fig. 1 is a schematic diagram of a measuring device according to a preferred 
embodiment; 

■ 

Fig. 2a and 2b show cross sections of a sample and two fiber tips in the situation 

wherein the mean free path of the photons is much larger than the diameter of the 
5 fibers; 

Fig. 3 shows the results of Monte Carlo simulations for a homogeneous medium; 
Fig. 4 a differential backscatter signal normalized at zero absorption for several 
scattering coefficients; 

Fig. 5 shows a differential backscatter signal of a dilute suspension of 0.2 ^m 
10 polystyrene spheres 

Fig. 6 shows the total differential backscatter signal as a function of the reflection 
coefficient in the range of 10-100 mm'^; 

Fig. 7 shows the measured and calculated avemge path length r as a function of the 
average scattering coefficient; 
IS Fig. 8 is a graph of measurement results for three different absorption coefficients ^Xa 
showing the dc-fiber signal /, the c^fiber signal J and the differential backscatter signal 
Khs ^ ^ function of wavelength; 

Fig. 9 shows a typical spectrum of an absorption curve A along with the specific 
absorption coefficient of Evans Blue dye; 
20 Fig. 10 shows a measiured ^ * as a function of the absorption coefficient ^ at A = 600 
nm; 

Fig. 1 1 shows typical spectra measured in a suspension of 1.0 (xm polystyrene spheres 
with and without Evans Blue dye; 

Fig. 12 graphically shows a molar extinction coefficient as a function of wavelength; 
25 Fig. 13 shows in vivo measurements and a fit of the differential backscatter signal Rhs 
in a human trachea realized usmg a fiber diameter of 400 pm, and 
Fig. 14 shows in vivo measurements and a fit of the differential backscatter signal R^s 
in a human trachea showing very low oxygenation indicative for lung tumor. 

30 A schematic diagram of a preferred embodiment according to the invention is 

shown in figure 1. The setup consists of a set of optical fibers for the deliveiy and 
collection of light to and from a sample 1 under investigation. Light fiom a Ught source 
2, for example a Tungsten Halogen lamp (Avantes HL-2000-FHSA), is led through a 
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first ann 3 of a bifurcated optical fiber. The bifurcated optical fiber is at a distal end 4 
coupled to a first distal end of a delivery- and-collection fiber 5 (in the following 
referred to as dc-fiber 5) which is small enough to be fit through a working channel of a 
clmical endoscope, not shown. A second distal end of the delivery- and-coUection fiber 
5 5 contacts the sample 1 . Alongside the dc-fiber 5, a collection fiber 6 is arranged to 
collect light reflected by the sample 1. The collection fiber 6 (referred to as c-fiber 6) is 
connected to a slave channel of a dual-channel spectrometer 7, for example an Avantes 
SD2000. Preferably, the dc-fiber 5 is polished at a small angle to reduce specular 
reflections. » 

1 0 Light reflected back fi:om the sample 1 into the c-fiber 6 is led directly into the 

slave channel of the dual-channel spectrometer 7. A second arm 8 of the bifurcated 
fiber is connected to a master channel of the dual-chaimel spectrometer 7. Light 
reflected into the dc-fiber 5 is coupled back into the bifurcated fiber, and reaches the 
dual-channel spectrometer 7 via the second aim 8 of the bifurcated fiber. An output of 

15 the spectrometer 7 is connected to an input of a processor 9 which is arranged to 
analyze signals &om the spectrometer 7. 

If only the dc-fiber 5 would be used to deliver and collect light to and from the 
sample 1, a large fraction of collected light is due to single backscattering fi:om small 
sample depths, see [1]. A single to multiple scattering ratio depends on the scattering 

20 coefficient and phase function of the sample 1 and on a diameter of the dc-fiber 5. The 
contribution of multiply scattered light to the signal of the dc-fiber 5 can be 
approxunately determined by combining the signal of the dc-fiber 5 with a signal 
coming from an additional fiber, i.e. the c-fiber 6 mentioned above. 

In [4] a differential backscatter signal Rts as a fimction of the wavelengdi X is 

25 determined using a formula like 



(1) 



30 



where 1(2.) is the signal from the dc-fiber 5 in contact wifli the sample 1, I„(X) is the 
signal bom the dc-fiber 5 submersed in a fluid with an appropriate refiactive index (for 
tissue: water would be appropriate), J^ute G) is the signal firom the dc-fiber 5 with the 
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probe-tip at a specific distance £ix>m a diffuse reflecting reference material with a large, 
preferably wavelength-independent reflectance coefficient (white speclialon) and/ji^ 
(k) is the signal from the dc-fiber 5 with the probe-tip at that same specific distance 
from a diffuse reflecting reference material with a small, preferably wavelength- 
5 independent reflectance coefficient (black spectialon). Furthermore, J (X) is tiie signal 
from the c-fiber 6 in contact with the sample 1 and JyMiack (X) is the signal from the c- 
fiber 6 with the probe-tip at the previously mentioned specific distance from the 
white/black spectralon. Fmally, c is a calibration constant fliat depends on the distance 
between tiie probe-tip and the reference materials. 
1 0 According to tiie invention, the processor 9 is arranged to calculate the physical 

feature using a predefined matiiematical model, tiie differential backscatter signal (Rts) 
and a curve fitting mechanism. In an embodiment, tiie diameter of tiie fibers 5, 6 are 
selected depending on a mean fi«e patii (m^) of photons sent into tiie sample 1. It is 
noted tiiat if tiie mean free patii can not be estimated before selecting a fiber diameter, 
15 mitiaUy two arbitiary fiber diameters may be selected. Afier curve fitting tiie measuring 
results using two different matiiematical models, it will show which model applies. 

Figures 2a and 2b depict fiber tips of tiie dc-fiber 5 and the c-fiber 6 in the 
situation wherem tiie mean fiee patii (m^) of photons coming out of tiie dc-fiber 5, is 
much larger tiian a diameter djiter of flie fibers 5, 6. In an embodunent, tiie diameters of 
20 botii fibers 5, 6 are of equal size, however it should be understood tiiat otiier selections 
are possible. In figure 2a, lines 21 and 22 show an example of a patii traveled by a 
detectable smgly scattered photon. In figure 2b. lines 23, 24, 25 and lines 23, 24, 26 
show two possible patiis of detectable multiply scattered photons. All multiple 
scattering events occur at such large distances from tiie fiber tip of tiie fibers 5, 6, tiiat 
25 tiie probability of detection of multiply scattered photons is roughly equal for tiie dc- 
fiber 5 and tiie c-fiber 6. The differential backscatter signal Ri„(X) will now purely be 
determined by singly scattered photons. 



In an embodiment, tiie diameter of tiie fibers 5, 6 are selected so tiiat m^><^ 
In tiie predefined matiiematical model of tiiis embodiment, tiie differential backscatter 
signal Rta(X) is an exponential fiinction of two times tiie mean free path. Below, an 
e2q}lanation for this model is given. 
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In the absQice of absorbers, the differential backscatter signal Rti(k) is 
proportional to the local, superficial scattering coefficient fh(^)=QscaW p ' As : 



— • ^Qi.i-pKyi,ii)'y^{yi)'p-As (2) 



where is an apparatus constant that depends amongst others on the distance 
between the probe tip and the reference materials (black and white spectralon), p(X, Q) 
is a function called the phase function where Q is the scattering angle, QscaW the 
scattering efficiency, p the concentration of substances present in the sample 1, and As 
10 the area of a scattering particle. For example, using a fused silica fiber with numerical 
aperture NA= 0.22, the differential backscatter signal Rbs(^) can be approximated by 

J 2ff 180 

RbsW^ Capp • — • Jd9 • Jde . sin(e) • piA,l 80) • //, (A) = Capp'- pQ, 1 iOy/j,(X) (3) 

0 170 



15 where ^ is the azdmuthal angle and 9 is the polar angle. 

Figure 3 shows the differential backscatter signal RtsCl) of measurements (see 
dots) of a dilute suspension of 0.2 urn polystyrene q)heres along with a calculation (see 
curve 32) according to Eq. (3). In figure 3 Ri,s(X) is shown using arbitrary units (a.u.). 

20 Also a value of Q„dar where Qradar = 4np(A,l SO)Q^(X), is indicated in the figure. 
Figure 3 shows excellent agreement between the measurement (i.e. the dots) and the 
calculation, which mdicates that if in^>rf^4er , the single scattering is indeed the 
donunant contributor to the differential backscatter signal Ri^) as defined in Eq. (1). 
A smgly scattered photon first travels fiom the tip of the dc-fiber 5 to a particle, 

25 and then (the same distance) back &om the particle to the tip of the dc-fiber 5 (or tip of 
c-fiber 6), see also figure 2a. So an average path length t(X) traveled by the measured 
single scattered photons is equal to two times the mean free path m^(X), i.e. 



30 



1(2.) =2-m£p(:i) 
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In the presence of n absorbing species with specific absorption coefficients fia^^^(X)y 
the differential backscatter signal becomes 

Rbsa) =Capp' •Pttl80)-//,(>l)-exp(- t(A)- 1; Pi • 
5 = Capp' • pCl 180) • mM) • exp(-2 • m^(>l)) • J) Pi • A.^ '^)) (5) 

where Capp' is an apparatus constant, p(A, 1 80) is the phase function, /^A) is the 
scattering coefficient of the medium, A is the wavelength of the first and second 
backscattered radiation, mfj^X) is the mean fi:ee palb as a function of the wavelength, n 

0 is the number of substances in the sample 1, pi is the concentration of absorber / present 
in a detection volume of the sample 1 , and /4z^^'(A) is the absorption coefficient of 
absorber / as a function of the wavelength. 

It is noted that in Eq. (S) the assumption is made that absorbers are 
homogeneously distributed and do not influence each other. The Eq. (S) may be 

S corrected for non-linear phenomena such as an inhomogeneous distribution of 
absorbers, see e.g. [8]. 

According to an embodiment, the specific absorption coefficients of the 
absorbers, the wavelength dependency of the scattering coefficient and the phase 
function /?, together with Eq. (5) are used in order to calculate the concentrations of all 

0 the absorbing substances present in the detection volume of the sample 1 . Since the 
detection volume is typically very small in the present invention, the extracted 
concentrations are highly spatially resolved. This is not possible with the known 
methods that are based on diffuse reflectance, and wherein the obtained concentrations 
are averages over large sample volumes, see e.g. [2]. 

5 The apparatus constant Capp' (Eq. 3) can be determined for a specific distance 

between the tip of the dc-fiber 5 and the reference materials (black and white 
spectralon). For a suspension of monodisperse polystyrene spheres of known size and 
concentration, the scattering coefficient /jg and the phase function p(l 80) can be 
calculated using Mie theory [4]. The apparatus constant Capp' simply follows from Eq. 

0 (3). In terms of the volume fraction / of the suspension, the radius of the spheres a and 
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the radar efficiency coefficient Gr<ttfar^)= 4n'pQ^ ,180)-Qjca(^) the apparatus constant is 
determined by 

Rb^X) « Capp'-p(X,180)-/^) = Capp*- 0.05968 • - • g^^) (6) 

a 

5 

According to another embodiment, the selected diameter dfiber is chosen so that 
the mean free path is smaller than dfiber- In this embodiment, the differential backscatter 
signal Rte is a function of the fiber diameter dfiber- This will be discussed in more detail 
below. 

10 When the mean free path of the photons is smaller than the selected fiber 

diameter (i.e. m^(A)<dfiber)9 the contribution of multiply scattered light to the 
dififerential backscatter signal RbsW of the single dc-fiber 5 cannot completely be 
removed using Eq. (1). In this case, it appears that the average path length of the 
photons contributing to the signal Rbs(^) becomes nearly independent of the optical 

15 properties of the sample 1. In this situation, multiple scattering events already occur at 
small distances from the tip of the dc-fiber 5. An analytical e2q}ression for the 
backscatter signal Rbs(^) is not available for this situation and Monte Carlo simulations 
were used to model the behaviour of Rbs(^) as a function of the diameter of the fibers 5, 
6 and of the optical properties of the sample 1. Figure 4 shows the results of Monte 

20 Carlo simulations using the MCML-code (Monte Carlo for Multi-Layered media) of 
Wang et al [6,7] for a homogeneous medium with an anisotropy value g=0.9. A flat 
circular incident beam with diameter dfii^r is directed onto the sample 1, and the 
differential backscatter signal is calculated by subtracting the total reflectance in the 
c-fiber 6 (with diameter dfiter and center located at a distance dfiber from the center of the 

25 incident beam) from the total reflectance in the dc-fiber 5 (with diameter dfiber 

overlapping the incident beam). Simulations were perfomied for sets of four different 
scattering coefficients (m=15, 25, 50 and 80 mm*^), four different fiber diameters 
(^6cr=200, 400, 600 and 800 jxm) and five different absorption coefficients (fia==0, 0.2, 
0.4, 0.6 and 0.8 mm"*). 

30 Figure 4 shows i^A* as a function of absorption coefficient pia where the open 

circles/dashed lines correspond to dfiifer=200 pm, the filled circles/dotted lines 
correspond to 4?&er=400 jmi, the open squares/solid lines correspond to dfiber=600 pm 
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and the filled squares/dash-dotted lines correspond to dfiber=SOO jun. The dififerential 
backscatter signal for each scattering coeflBcient iis was normalized to the 
(dfiber=200 |jm, jia=0 mm'^) case. Figure 4 shows that in the absence of absorption, Le. 
jOa = 0, the differential backscatter signal Rts depends linearly on the scattering 
5 coefficient Hs. Furthennore, the slope of the straight lines (signifying the relation 
between Rbs and jUg) depends only on the fiber diameter and is independent of the 
scattering coefficient pig. The latter is more clearly demonstrated in figure 5, where the 
differential backscatter signal Rts is normalized to unity at zero absorption for all 
scattering coefficients p^. The open circles correspond to ^^r=200 jjin and the filled 
10 squares correspond to dftber^SOO pm. These Monte Carlo simulations therefore suggest 
that in the situation where m^<dfiber, the diameter of the fibers 5, 6 determines the 
average path length x of the measured photons. The backscatter signal Ri,s for this range 
of parameters can thus be written as 

15 Rbs(Ji) = Cr m • exp(-T ■ pa) = Cr Hs * exp(-C2 • dfiter ' Pa) (7). 



where Ci and C2 are constants, 1; is the average path length, pa is the absorption 
coefficient, m is the scattering coefficient and d/iberis the fiber diameter of the fibers 5, 
6. 

20 An exact analytical ejqpression for Rt/X) is not available due to the large 

contribution of multiple scattering events to the signal. Measurranents were done for 
determining a total integrated backscatter signal Rtot for a range of A between 400-900 
nm, using the formula 

900inii 

25 RtotGi«)= J dXRbsd/j,) (8) 

Figure 6 shows that the integrated total backscatter signal Riot(jLis} is proportional to 
MsW in the relevant range of jUg between 10-100 mm"^ Therefore, in the absence of 
absorbers, it follows that 
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RbM = C^'/isG) (9) 

which is in agreement with the Monte Carlo simulations. 

In the presence of « absorbing substances in a suspension, with specific 
absorption coefficients fi^^'^'iX), the dijBferential backscatter signal becomes 



= -expC-r • 2] Pi ' tJ^'^m (10) 



where t is the average path length of the detected backscattered photons and pi is the 
concentration of the substance /. 

Non-linear phenomena such as an inhomogeneous distribution of absorbers are 
not incorporated in Eq. (10), but can be added by the skilled person, see e.g. [8]. 

1 5 Figure 7 shows the measured and calculated average path length r as a function 

of the average scattering coefficient <HsQ)> (with 500 nm < A, < 700 mn) for dfiUr = 0.4 
mm and for absorption coefficient Ha(X) = 2.0 mm"* at ;i=600 mn. In figure 7 
measurement are indicated by dots, the r = 2-m^ curve is indicated by a line 71 and 
Monte Carlo simulations are depicted by dashed lines. Identical results were obtained 

20 for suspensions with an absorption coefficient of /4, = 1 .0 mm"^ at 600 nm. The average 
path length r was determined using suspensions of polystyrene spheres with different 
sizes and concentrations to vary the scattering coefficient Hs(X). The anisotropy g of 
these suspensions was in the range of 0.8-0.9. Evans Blue dye was added as an 
absorber, and the average path length r was calculated fiom Eqs. (9) and (10) and 

25 knowledge of the concentrations and specific absorption coefficient of Evans Blue, as 
will be known to the skilled person. 

Looking at the measured average path lengths of figure 7, it clearly shows that for 
large scattering coefficients {jjk = 10-100 mm'\ the range relevant for tissue) the 
average path length r is independent of the scattering coefficient /4 to within 10% and 

30 approxmiately equal to half the fiber diameter (r « 0.24 mm while dfih^ = 0.40 mm), in 
agreement with the Monte Carlo simulations (the dashed lines correspond to Monte 
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Carlo calculations for djiter = 0.2, 0.4, 0.6 and 0.8 mm). For small scattering coefficients 
(e.g. A < 5 mm*^) the average path length r is well described by r= 2-m§>, according 
to Eq.(4), see line 71 . Figure 7 also clearly demonstrates that the transition from the 
•single scattering regime' to the 'constant path length regime' occurs for mean free 
5 paths of the order of the fiber diameter. It is therefore expected that single scattering 
will prevail over a larger range of scattering coefficients for fiber diameters smaller 
than 400 \mx. 

In the following, the effect of absorption on the average path lei^ r will been 
examined in more detail. Various concentrations of Evans Blue dye were added to a 

10 suspension of polystyrene spheres with scattering coefficients of 35 mm"^ The 

concentrations of Evans Blue (EB) dye were varied such that the absorption coefficient 

at 600 nm was in the range of 0 to 2 mm"^ Typical results of the differential 
backscatter signal Rbs for three different absorption coefficients /Ja are shown in figure 
8. Note that the signal I of the dc-fiber 5 is plotted on a different vertical scale than the 

1 5 signal J of the c-fiber 6 and the differential backscatter signal Rbs- 

The spectra with Evans Blue present in the suspension were divided by the 
spectrum with no Evans Blue present and the negative nattiral logarithm of the ratio 
R^^/R^ was determined: 

20 ^ = -hi(R^^/R°) = T-p-/4,^'^ (11) 



where p is the concentration of Evans Blue, and ^4,^^^^'^^ is the specific absorption 
coefficient of Evans Blue. 

25 Figure 9 shows a typical spectrum of an absorption curve 92, along with the 

specific absorption coefficient fj^"^ of Evans Blue dye, see curve 94. 

For all concentrations, an area^i* under the absorption curve 92 was determined 
in the wavelength range X between 500 and 650 nm. From Eq. (1 1) it follows that if the 
average path length x is independent of the absorption coefficient /Ja^"", area ^ * should 
30 depend linearly on the concentmtion p of the species m the suspension. 

Figure 10 shows a measured ^ * as a function of the absorption coefficient /4, at 
600 nm. A curve fitting is done, resulting in a line 104 for fds(X )=35 mm"^ Figure 10 
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shows that the average padi length t is indeed independent of the absorption coefficient 
/ib in the range 0-2 mm"^ . 

From the previous results of figure 4 to 10, it diows that for ni^< dfihen the 
dijfferential backscatter signal Rts is described by Eq. (7) with C2 » 0.6. Figure 1 1 
shows typical spectra measured in a suspension of 1 .0 |Jiu polystyrene spheres with and 
without Evans Blue dye (|ia= 2 and 0 mna'^ at 600 nm, respectively). From Eq. (7) and 
figure 7, the relation between the differential backscatter signals with and without 
absorber is given by 



10 Rbs(^,t^ = Rts(^0yexp('0,24^iJ^ (12) 

The calculated spectrum according to Eq. (12) is plotted as a dashed line 1 10 in 
figure 1 1 and shows excellent agreement with the measured Rt^A^^j), see line 1 1 1 . In 
figure 1 1 , line 1 12 depicts Rb^^ 0). 
15 In short, the average path length t of photons measured \^en subtracting the 

signals of the c-fiber 6 from the dc-fiber 5 using Eq. (1), is independent of the optical 
properties of the sample 1 and approximately equal to half the diameter of the fibers 5, 
6 used, as long as the fiber diameter is larger than the m^. 

20 In a specific embodiment, the device according to the invention is arranged to 

determine concentrations of oxygenated blood in tissue. Since the scattering coefficient 
of tissue /^'"^ is in the range of 10-100 mm'^ the fiber diameter should be smaller 
than a certain maximum diameter dmax where d^ax is between 10 and 100 |im, so for 
example smaller than 50 imi, in order to measure predominantly single scattering in 

25 tissue. In this case, Eq. (5) holds. For fibers 5, 6 with much larger diameters (e.g. 200 
or 400 |im), the differential backscatter signal Rb/X) is described by Eq. (10) with ra* 

O.&dfiber^ 

It is presently known that the wavelength dependence of the scattering coefficient 
in tissue ^^'^ can be adequately described by an empirical power-law fimction, see 
30 also [3], [4], [5]. 

^iissueg^y. a • (13) 
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with a and b constants that depend on the size, concentration and relative refractive 
index of the scatterers (Le. substances) present in the detection volume. 
The domioant absorbers in tissue in the visible wavelength range are oxygenated and 
deoxygenated blood. Thus in tissue Eq. (10) becomes 

Rbl^) = Capp • aA,-^ • exp(. Q.&dfiber • Pblood • (So2 ' Ma'''^'^ + (I-S02) ' /^^^^ )) 

= C'app • A,^ • exp(- 0.6'rf^er ' Pblood * (Soz ' /ii,^^'^ + (I-S02) " Pa^''''''' )) (14) 



where pwood is flie concentration of blood, S02 is the blood oxygenation (percentage 
oxygen saturation) in a certain detection volume, Capp is a constant that depends on the 
caUbration constant c, Capp is C^p-a, A. is the wavelength, b is the slope of the 
scattering coef&cient defined in Eq. 13, fii^^^^ is the specific absorption coefficients of 

15 fiilly oxygenated blood, jj^^'^^ is the specific absorption coefficients of fully 
deoxygenated blood. 

Non-linear phenomena such as an inhomogeneous distribution of absorbers are 
not incorporated in Eq. (14), but can be added by the skilled person, see e.g. [8], 

Since the specific absorption coefficients of fully oxygenated (pj^^^^ and fully 

20 deoxygenated {jjj^^*^^^ blood are well known, see figure 12, Eq. (14) can be fitted to 
the measured data to yield the slope b of the scattering coefficient i4^^^ the 
concentration pbiood and the oxygen saturation S02 of the blood present in the detection 
volume. When a correction is made for the inhomogeneous distribution of blood in the 
vessels, a vessel diameter D may be determined as well. Since the average detection 

25 depth is small (e.g. 0. 1 mm), the blood present in the detection volume when measuring 
non-invasively is located in capillaries. 

In figure 13, in vivo measurements of backscattering in a human trachea together 
with a fit using Eq. (14) are shown. The measurements are realized using a fiber 
diameter of 400 pm. The dots depict the measurements and curve 130 is a fitting curve. 

30 In figure 13,6 = -0.94 and the oxygenation Sq2 = 95%. 

The present invention can be used for tumor detection. Tumor growth may, due 
to its excessive oxygen consumption, be accompanied by a low capillary oxygen 
saturation, which can only be assessed using a very localized measurement. Since (pre- 
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)cancerous tissue is generally more heterogeneous than normal tissue, the standard 
deviation of multiple measurements is likely to be larger for (pre-)cancerous tissue than 
for normal tissue. Standard deviations in the measurements can be calculated for the 
oxygen saturation, the blood concentration, the blood vessel diameter, and the slope b 
5 of the scattering coeflBcient juj^. It is noted that the invention is by no means 
restricted to determine a concentmtion of a substance as the physical feature. All 
features, mentioned in the previous phrase can be regarded as physical features. 

An example of a measiirement of a lung tumor is shown in figure 14. The shape 
10 of the dip in the wavelength range of 500-600 nm in this figure demonstrates the 
depletion of oxygen from the capillaries of this tumor due to its excessive oxygen 
consumption. 

When a needle-probe is used, the local oxygenation and scattering coefficient 
juj'^^ can be measured invasively . This could be helpfiil in determining tumor-margins 

1 5 intra-operatively in real-time, for instance during resection of a breast-tumor. 

According to an embodiment, the device comprises multiple probes and a 
multichaimel spectrometer for multiple simultaneous measurements on dififerent 
locations of the sample 1. Using this device, multiple measurements can be made 
simultaneously on different locations of for example a suspicious lesion. 

20 In yet another embodiment, the device comprises at least two pairs of fibers, 

having different fiber diameters. For example, when a pair of fibers with 100 |im, a pair 
of fibers with a diameter of 200 pm and a pair of fibers with a diameter of 400 pm are 
used, information from different depths in the sample 1 can be obtained as the avemge 
path length increases with increasing fiber diameter. 

25 The method and appamtus according to the invention can also be used to analyze 

local drug concentmtions. From Eq. (10) it follows that if the specific absorption 
coefficient of a certain drug is known, the local concentration p of that substance can 
be determined using the invention. 

Another possibility of the present invention is to monitor glucose concentrations. 

30 The scattering coefficient ///^ depends among others on the relative refi:active index 
of the scatterers with respect to the surrounding medium (in tissue: cytoplasm). The 
refi:active mdex of the surrounding cytoplasm is likely to depend on the concentration 
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of glucose. A change in the glucose concentration will therefore likely affect the slope 
b of the scattering coefiScient ^^^^ , see Eq. (13). 

■ 

Whilst specific embodiments of the invention have been described above, it will 
be appreciated that the invention may be practiced otherwise than as described. For 
5 example, a concentration of a substance in polluted water may be calciilated. The 
description is not intended to limit the scope of the invention. 
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